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Abstract. Radiative decays of the B meson are known to provide important constraints on the 
MSSM and many other realistic new physics models in the sub-TeV range. The inclusive branching 
ratio 33(B — > X s y) being the key observable is currently measured to about ±7% accuracy. Reaching 
a better precision on the theory side is a challenge both for the perturbative QCD calculations and 
for analyses of non-perturbative hadronic effects. The current situation is briefly summarized here. 1 



The loop-generated effective bsy coupling is sensitive to new physics at scales 
(a few)x ^(lOOGeV) even in the simplest theories with Minimal Flavor Violation, 
like the Two-Higgs-Doublet Model or the Minimial Supersymmetric Standard Model. 
Measuring the inclusive 38 (B — > X s y) is the most efficient way to constrain this coupling. 
The Standard Model (SM) contribution to this branching ratio forms a background to 
effects we would like to put bounds on. Its precise QCD calculation has been a chal- 
lenge for a long time. A satisfactory situation has not yet been reached even for the 
perturbative contributions. 

The current experimental world averages (for E r > 1.6GeV in the decaying meson 
rest frame) read: 



They have been obtained by combining the measurements of CLEO [3], BABAR [4]-[6] 
and BELLE [7, 8] with different lower cuts Eq on the photon energy, ranging from 1.7 
to 2.0 GeV. An extrapolation in Eq down to 1 .6 GeV has been performed simultaneously. 
Uncertainties due to modeling the photon energy spectrum that matter both for the 
averaging and extrapolation have been singled out in Eq. (1). Ref. [1] gives a larger error 
than [2] because it uses results at Eq > 1.8 GeV from the older measurements [3]-[7] 
only, not taking into account the most precise ones from Ref. [8]. 

Calculations including &(aj) and ff(a^ m ) effects in the SM give [9, 10] 
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INTRODUCTION 



38{B X s y) 




(3.55±0.24 exp ±0.09 model ) x 10~ 4 [1], 
(3.50±0.14 exp ±0.10 model ) x 10- 4 [2]. 



(1) 



38(B^X s y) = (3.15 ±0.23) x 10 
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(2) 
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where the error is found by adding in quadrature the non-perturbative (5%), perturbative 
(3% + 3%) and parametric (3%) uncertainties. The result in Eq. (2) is consistent with the 
averages (1) at the 1.2a level. Its evaluation is based on an approximate equality of the 
hadronic and perturbatively calculable partonic decay widths 

T(B X s y) Er>Eo ~ T(b XP r ) Ey>Eo , (3) 

where X[ stands for s, sg, sgg, sqq, etc. This approximation works well only in a certain 
range of Eq, namely when Eq is large {Eq ~ mi,/ 2) but not too close to the endpoint 
(m b — 2Eq 3> A QCD ). Corrections to Eq. (3) of various origin have been widely discussed 
in the literature, most recently in Ref. [11]. 



PERTURBATIVE CALCULATIONS 

Radiative B decays are most conveniently analyzed in the framework of an effective 
theory that is obtained from the SM by decoupling of the W boson and all the heavier 
particles. Flavor-changing weak interactions are then given by 

^weak~ Y,Q(H)Q U (4) 

i 

where the operators Qj are built of the light fields only, while Q{fl) are the Wilson coeffi- 
cients. Once the higher-order electroweak and/or CKM-suppressed effects are neglected, 
eight operators matter for B — > X s y. They are displayed in Fig. 1. Their Wilson coeffi- 
cients at the scale fit, ~ m^/2 are presently known up to the Next-to-Next-to-Leading 

Order (NNLO) in QCD, i.e. including corrections up to [a] [a s \a ^ j j 

The necessary matching [12, 13] and anomalous dimension [14, 15, 16] calculations 
involved Feynman diagrams up to three and four loops, respectively. 

Once the Wilson coefficients are at hand, the partonic decay rate is evaluated accord- 
ing to the formula 

8 

r(b^X?y) Er>Eo =N£ Q(Hb)Cj {Hb)Gij (Eq, fi b ) , (5) 
i'J=i 

where N = \V*V tb \ 2 {G 2 F m 5 b a &m ) / {32n A ) . At the Leading Order (LO), we have Gy = 
8a dp, while the (?{a s ) Next-to-Leading Order (NLO) contributions are known since a 




FIGURE 2. Examples of Feynman diagrams that contribute to G77, G78 and G21 at 0(0%). Dashed 
vertical lines mark the unitarity cuts. 



long time (see Ref. [17] for a description and references). 2 At the NNLO, it is sufficient 
to restrict our attention to z',j e {1,2,7,8} because the penguin operators have very 
small Wilson coefficients ({Cs^fii^l < |Q(^)| ~ a s (jlj,)/K). In the following, we 
shall treat the two similar operators Q\ and Q2 as a single one (represented by Q2), and 
consider six independent cases of the NNLO contributions to G,y. 

Three of those six cases (G77, G78 and G27) involve the photonic dipole operator 
Qj. Examples of the corresponding contributions to the decay rate are shown in the 
subsequent columns of Fig. 2 as propagator diagrams with unitarity cuts. Two-, three- 
and four-body final states appear in the first, second and third rows, respectively. The 
rows cannot be considered separately because cancellation of IR divergences takes place 
among them. While G77 was found already several years ago [19]— [23], the complete 
calculation of G78 has been finalized only very recently [24, 25]. Evaluation of G27 is 
still in progress (see below). 

The remaining three cases (G22, G28 and G&g) receive contributions from diagrams 
like those displayed in Fig. 3. Diagrams in the first row involve two-body final states 
and are IR-convergent. They are just products of the known NLO amplitudes. Three- 
and four-body final state contributions remain unknown at the NNLO beyond the 
BLM approximation [26]. The BLM calculation for them has been completed very re- 
cently [27, 28] providing new results for G28 and Ggg, and confirming the old ones [29] 
for G22- The overall NLO + (BLM-NNLO) contribution to the decay rate from three- 
and four-body final states in G22, G28 and Ggg remains below 4% due to the phase-space 
suppression by the relatively high photon energy cut Eq. Thus, the unknown non-BLM 
effects here can hardly cause uncertainties that could be comparable to higher-order 
0(cts) uncertainties in the dominant terms (G77 and G27). One may conclude that the 
considered G ( - y - are known sufficiently well. 

It follows that the only contribution that is numerically relevant but yet unknown at 



These statements at the LO and NLO hold up to tiny but yet unknown contributions involving four-quark 
penguin operators and syqq final states. See Appendix E in Ref. [18]. 




FIGURE 3. Examples of Feynman diagrams that contribute to G22, (?28 and Gss at 0(0%). 



the NNLO is G27. So far, it has been evaluated for arbitrary m c in the BLM approxi- 
mation [30, 29] supplemented by quark mass effects in loops on the gluon lines [31]. 
Non-BLM terms have been calculated only in the m c 3> m^/2 limit [10, 32], and then 
interpolated downwards in m c using BLM-based assumptions at m c = 0. Such a pro- 
cedure introduces a non-negligible additional uncertainty to the calculation, which has 
been estimated at the level ±3% in the decay rate. 

As a first attempt to improve the situation, a calculation of G27 at m c = has 
been undertaken [33, 34]. Two- and three-particle cut contributions have been already 
found [35]. They contain an IR divergence which should be canceled by diagrams with 
four-particle cuts [34] . 

A recently started calculation [36] for arbitrary m c is supposed to cross-check the 
m c = result and, at the same time, make it redundant, because no interpolation in 
m c will be necessary any more. The method to be used is the same as in the BLM 
calculation of Ref. [31]. However, the number of master integrals to be considered is 
now much larger (^(500)). A system of differential equations for them with respect 
to the variable z = ml/ml needs to be solved (numerically) along an ellipse in the 
complex z-plane. The boundary conditions at z ^> 1 are going to be found with the help 
of asymptotic expansions. The most computer-power demanding part is the integration- 
by-parts reduction to master integrals that is currently being performed. 

NON-PERTURBATIVE CONTRIBUTIONS 

The question to what accuracy the approximate equality (3) holds has been subject 
of many investigations since early 1990's. However, a quantitative analysis of all the 
dominant contributions to the resulting uncertainty in £${B — > X s y) has been performed 
only very recently [11]. 

Obviously, corrections to Eq. (3) depend on Eq. As already mentioned in the intro- 
duction, they are minimized at a certain "optimal" value of Eq that is high enough 
(Eq ~ m^/2) but not too close to the endpoint (mt, — 2£b 3> A). The value of Eq = 
1.6 GeV ~ m h /3 that has been chosen in Ref. [18] as default seems to have a chance to 




FIGURE 4. Examples of diagrams describing non-perturbative contributions to the "27" interference 
term due to soft gluons originating from the B-meson initial state. 



be in the vicinity of the optimal point. In the following, I will discuss non-perturbative 
effects at this value of the cutoff, leaving aside the problem of photon energy extrapola- 
tion in the experimental averages. 3 

So long as only the photonic dipole operator Qj is considered, non-perturbative 
corrections to Eq. (3) for — 2Eq 3> A can be described in terms of the so-called fixed- 
order approach that has been derived [37] using the optical theorem and the Operator 
Product Expansion. The corrections can then be written as a series in (A/m^)"af with 
n = 2, 3,4, . . . and k = 0, 1,2, . . ., where perturbatively calculable coefficients multiply 
matrix elements of local operators between the 5-meson states at rest. Such matrix 
elements (at least the leading ones) can be extracted from measurements of observables 
that are insensitive to new physics, like the semileptonic B — > X c ev decay spectra or 
mass differences between various ^-flavored hadrons. Coefficients at the terms of order 
A 2 /ml and A 3 /m| have been evaluated in Refs. [38, 39] and [40], respectively. Very 
recently, a calculation at order a s A 2 /ml has been completed [41]. The result explodes 
near the endpoint Eq ~ but remains perfectly consistent with the fixed-order 

approach at Eq = 1.6 GeV or even 1.7 GeV. Thus, non-perturbative corrections to the 
"77" interference term are well under control. In many phenomenological analyses, 
normalization to the semileptonic rate is applied in such a way that most of those 
corrections cancel, leaving out a sub-percent effect. 

The most important non-perturbative uncertainty originates from the "27" interference 
term (that stands for "27" and "17"). In Ref. [11], photons that can be treated in analogy 
to the "77" term are called "direct", while all the other ones are called "resolved", i.e. 
produced far away from the b-quark annihilation vertex. Contributions from the resolved 
photons can still be written in terms of a series in powers of (A/m^)"af , but this time the 
(n = 1 , k = 0) term is non- vanishing when m c is treated as & ( V / Am^) . Moreover, they are 
uncertain, as they depend on matrix elements of non-local operators that cannot be easily 
extracted from other measurements. Diagrams representing such terms are displayed in 
Fig. 4, where the external gluon is understood to be soft, while the other one (if present) 
is considered to be non-soft. 

If the charm quark was heavy enough (m^/m^ ^> A), its loop in the first diagram of 
Fig. 4 would become effectively local for soft gluons, and we would be back to the local 
operator description, as in the "77" term. This limit has been analyzed in Refs. [42]-[46] . 



J Measurements at Eq = 1.6 GeV or 1.7 GeV will always be plagued with much larger background 
subtraction errors than the ones at 1 .9 GeV or so. Thus, we will always need to search for a proper balance 
between those errors and uncertainties due to the photon spectrum modeling. Background subtraction 
requires modeling, too. 



A series of the form 



n=0 



was found as a relative correction to Eq. (3). Explicit results for all the coefficients b n 
showed that they are small and quickly decreasing with n, which led to a conclusion 
that the first term in the series is a good approximation to the whole correction even in 
the m^A/m^ ~ ^(1) case that we encounter in Nature. The leading term is proportional 
to a local operator matrix element that can be extracted from the measured B-B* mass 
difference. This way, a relative correction of around +3% to Eq. (3) has been found. 

This conclusion has recently been questioned in Ref. [11] on the basis of realistic 
shape function models that allowed to vary m c in the physically interesting range, and 
test applicability of the expansion (6). It has been found that the first term of such an 
expansion in not really a good approximation if we allow for alternating-sign subleading 
shape functions (see Eq. (108) in that paper). Alternating signs were necessary to 
overpass normalization constraints, and make the shape function exponential tails wide 
enough for more energetic (though still soft) gluons. This is the main source of the 
overall ±5% non-perturbative uncertainty in the branching ratio that was estimated 
in Ref. [11]. So long as m c is treated as 0(y/Amb), the considered correction is just 
&(A/mb). Other (smaller) corrections studied in that paper were of order &(a s A/mb), 
and did not originate from the "27" interference term. 

Apparently, the 0{a s A./mb) corrections alone were the reason for assigning a ±5% 
non-perturbative uncertainty to the branching ratio in Refs. [9, 10]. Thus the two ±5% 
uncertainty estimates agree just by coincidence. The main worry in Refs. [9, 10] were 
the second and third diagrams in Fig. 4. Although they look like a s corrections to the 
first one, they are not necessarily unimportant given that the smallness of the first one is 
partly accidental. It would be interesting to test their relevance using the shape function 
methods. This would make the analysis of &(a s K/mb) uncertainties really complete. 

In the end, let us recall that there exist non-perturbative corrections to Eq. (3) that are 
not suppressed by A/ra^ at all. Their intuitive description can be found in Ref. [47]. In 
particular, collinear photon emission effects belong to this class [48, 27]. Fortunately, 
they are numerically small due to interplay of several minor suppression factors. 



SUMMARY 

Given the present consistency of measurements and SM calculations, observing clean 
signals of new physics in B — > X s y is unlikely, even if the uncertainties are reduced by 
factors of 2 on both sides, which may be hoped for in the Super-5 era. However, achiev- 
ing such a reduction is worth an effort, as it would lead to strengthening constraints on 
most popular beyond-SM theories. Several new perturbative NNLO results have been 
published this year, and new ones are expected in 2011. As far as the non-perturbative 
corrections are concerned, they are still dominated by unknown contributions, but at 
least their estimates are now based on calculations rather than order-of-magnitude con- 
siderations. Some room for improvement seems to remain both in the &(h/m,b) and 
&{a s A./mb) cases. 



ACKNOWLEDGMENTS 



This work has been supported in part by the Ministry of Science and Higher Education 
(Poland) as research project N N202 006334 (2008-11) and by the EU-RTN program 
"FLAVIAnet" (MRTN-CT-2006-035482). Partial support from the DFG through the 
"Mercator" guest professorship program is gratefully acknowledged. 



REFERENCES 



1. D. Asner et al. (Heavy Flavor Averaging Group), arXiv: 1010. 1589. 

2. M. Artuso, E. Barberio and S. Stone, PMC Phys. A 3, 3 (2009) [arXiv:0902.3743]. 

3. S. Chen et al. (CLEO Collaboration), Phys. Rev. Lett. 87, 251807 (2001) [hep-ex/0 108032]. 

4. B. Aubert et al. (BABAR Collaboration), Phys. Rev. D 72, 052004 (2005) [hep-ex/0508004]. 

5. B. Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 97, 171803 (2006) [hep-ex/060707 1 ] . 

6. B. Aubert et al. (BABAR Collaboration), Phys. Rev. D 77, 051 103 (2008) [arXiv:071 1.4889]. 

7. K. Abe et al. (BELLE Collaboration), Phys. Lett. B 511, 151-158 (2001) [hep-ex/0 103042]. 

8. A. Limosani et al. (Belle Collaboration), Phys. Rev. Lett. 103, 241801 (2009) [arXiv:0907.1384]. 

9. M. Misiak et al, Phys. Rev. Lett. 98, 022002 (2007) [hep-ph/0609232]. 

10. M. Misiak and M. Steinhauser, Nucl. Phys. B 764, 62-82 (2007) [hep-ph/0609241]. 

11. M. Benzke, S. J. Lee, M. Neubertand G. Paz, JHEP 1008, 099 (2010) [arXiv: 1003.5012]. 

12. C. Bobeth, M. Misiak and J. Urban, Nucl. Phys. B 574, 291-330 (2000) [hep-ph/99 10220]. 

13. M. Misiak and M. Steinhauser, Nucl. Phys. B 683, 277-305 (2004) [hep-ph/0401041]. 

14. M. Gorbahn and U. Haisch, Nucl. Phys. B 713, 291-332 (2005) [hep-ph/041 1071]. 

15. M. Gorbahn, U. Haisch and M. Misiak, Phys. Rev. Lett. 95, 102004 (2005) [hep-ph/0504194]. 

16. M. Czakon, U. Haisch and M. Misiak, JHEP 0703, 008 (2007) [hep-ph/06 12329]. 

17. A. J. Buras and M. Misiak, Acta Phys. Polon. B 33, 2597-261 1 (2002) [hep-ph/0207 1 3 1 ] . 

18. P. Gambino and M. Misiak, Nucl. Phys. B 611, 338-366 (2001) [hep-ph/0 104034]. 

19. K. Melnikov and A. Mitov, Phys. Lett. B 620, 69-79 (2005) [hep-ph/0505097]. 

20. I. R. Blokland, A. Czarnecki, M. Misiak, M. Slusarczyk and F. Tkachov, Phys. Rev. D 72, 033014 
(2005) [hep-ph/0506055]. 

21. H. M. Asatrian, A. Hovhannisyan, V. Poghosyan, T. Ewerth, C. Greub and T. Hurth, Nucl. Phys. B 
749, 325-337 (2006) [hep-ph/0605009]. 

22. H. M. Asatrian, T. Ewerth, A. Ferroglia, P. Gambino and C. Greub, Nucl. Phys. B 762, 212-228 
(2007) [hep-ph/0607316]. 

23. H. M. Asatrian, T. Ewerth, H. Gabrielyan and C. Greub, Phys. Lett. B 647, 173-178 (2007) [hep- 
ph/0611123]. 

24. H. M. Asatrian, T. Ewerth, A. Ferroglia, C. Greub and G. Ossola, arXiv: 1005.5587. 

25. T. Ewerth, Phys. Lett. B 669, 167-172 (2008) [arXiv:0805.391 1]. 

26. S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 28, 228-234 (1983). 

27. A. Ferroglia and U. Haisch, arXiv: 1009.2144. 

28. M. Misiak and M. Poradzihski, arXiv: 1009.5685. 

29. Z. Ligeti, M.E. Luke, A.V. Manohar and M.B. Wise, Phys. Rev. D 60, 034019 (1999) [hep- 
ph/9903305]. 

30. K. Bieri, C. Greub and M. Steinhauser, Phys. Rev. D 67, 1 14019 (2003) [hep-ph/0302051]. 

31. R. Boughezal, M. Czakon and T. Schutzmeier, JHEP 0709, 072 (2007) [arXiv:0707.3090]. 

32. M. Misiak and M. Steinhauser, Nucl. Phys. B 840, 271-283 (2010) [arXiv: 1005.1 173]. 

33. R. Boughezal, M. Czakon and T. Schutzmeier, in preparation. 

34. M. Czakon, P. Fiedler, T. Huber and T. Schutzmeier, in progress. 

35. T. Schutzmeier, Matrix elements for the B — > X s y decay at NNLO, Ph.D. Thesis, University of 
Wiirzburg, 2009. 

36. M. Czakon, R. N. Lee, M. Misiak, A. V. Smirnov, V. A. Smirnov and M. Steinhauser, in progress. 

37. J. Chay, H. Georgi and B. Grinstein, Phys. Lett. B 247, 399-405 (1990). 



38. I.I. Bigi et al, "A QCD 'manifesto' on inclusive decays of beauty and charm," DPF92 Proceedings 
(Batavia, November 1 992) [hep-ph/92 1 2227] . 

39. A. F. Falk, M. E. Luke and M. J. Savage, Phys. Rev. D 49, 3367-3378 (1994) [hep-ph/9308288]. 

40. C. W. Bauer, Phys. Rev. D 57, 5611-5619 (1998) [Erratum-ibid. D 60, 099907 (1999)] [hep- 
ph/9710513]. 

41. T. Ewerth, P. Gambino and S. Nandi, Nucl. Phys. B 830, 278-290 (2010) [arXiv:091 1.2175]. 

42. G. Buchalla, G. Isidori and S.J. Rey, Nucl. Phys. B 511, 594-610 (1998) [hep-ph/9705253]. 

43. Z. Ligeti, L. Randall and M.B. Wise, Phys. Lett. B 402, 178-182 (1997) [hep-ph/9702322]. 

44. A.K. Grant, A.G Morgan, S. Nussinov and R.D. Peccei, Phys. Rev. D 56, 3151-3154 (1997) [hep- 
ph/9702380]. 

45. M.B. Voloshin, Phys. Lett. B 397, 275-278 (1997) [hep-ph/96 12483]. 

46. A. Khodjamirian, R. Riickl, G. Stoll and D. Wyler, Phys. Lett. B 402, 167-177 (1997) [hep- 
ph/9702318]. 

47. M. Misiak, Acta Phys. Polon. B 40, 2987-2996 (2009) [arXiv:0911.1651]. 

48. A. Kapustin, Z. Ligeti and H.D. Politzer, Phys. Lett. B 357, 653-658 (1995) [hep-ph/9507248]. 



